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Abstract 

A few years ago (see [1]) two of us introduced, motivated by the 
study of certain forest-fire processes, the self-destructive percolation 
model (abbreviated as sdp model). A typical configuration for the 
sdp model with parameters p and S is generated in three steps: First 
we generate a typical configuration for the ordinary percolation model 
with parameter p. Next, we make all sites in the infinite occupied 
cluster vacant. Finally, each site that was already vacant in the be- 
ginning or made vacant by the above action, becomes occupied with 
probability 5 (independent of the other sites). 

Let d{p, S) be the probability that some specified vertex belongs, 
in the final configuration, to an infinite occupied cluster. In our ear- 
lier paper we stated the conjecture that, for the square lattice and 
other planar lattices, the function 6{-, •) has a discontinuity at points 
of the form {pc,d), with 6 sufficiently small. We also showed (see [2]) 
remarkable consequences for the forest-fire models. 

The conjecture naturally raises the question whether the function 
9{-, •) is continuous outside some region of the above mentioned form. 
We prove that this is indeed the case. An important ingredient in our 
proof is a (somewhat stronger form of a) recent ingenious RSW-like 
percolation result of BoUobas and Riordan ([4]). 

1 Introduction and outline of results 
1.1 Background and motivation 

The self-destructive percolation model on the square lattice is described as 
follows: First we perform independent site percolation on this lattice: we 
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declare each site occupied with probabihty p, and vacant with probability 
I — p, independent of the other sites. We will use the notation {V ^ W} 
for the event that there is an occupied path from the set of sites V to the 
set of sites W. We write {V 00} for the event that there is an infinite 
occupied path starting at V. 

Let, as usual, 9{p) denote the probability that a given site, say O = (0,0), 
belongs to an infinite occupied cluster. It is known that there is a criti- 
cal value < Pc < 1 such that 9{p) > for all p > Pc, and 6(j)) = for 
all p < Pc- Now suppose that, by some catastrophe, the infinite occupied 
cluster (if present) is destroyed; that is, each site in this cluster becomes 
vacant. Further suppose that after this catastrophe we give the sites inde- 
pendent 'enhancements', as follows: Each site that was already vacant in 
the beginning, or was made vacant by the catastrophe, becomes occupied 
with probability 5, independent of the others. Let Vp^s be the distribution 
of the final configuration. 

A more formal, and often very convenient description of the model is as 
follows: Let Xi, i e I? he independent — 1 valued random variables, each 
Xi being 1 with probability p and with probability 1 — p- Further, let Yi, 
i € 1?., be independent 0—1 valued random variables, each Yi being 1 with 
probability 6 and with probability 1 — 5. Moreover, we take the collection 
of YiS independent of that of the Xj's. Let X*,i G 1? be defined by 




1 if Xj = 1 and there is no X- occupied path from z to 00 
otherwise, 



(1) 

where by 'X-occupied path' we mean a path on which each site j has Xj = 1. 
Finally, define Zi = X* V Y^. This collection {Zi,i £ I?) is (with meaning 
'vacant' and 1 'occupied') what we called 'the final configuration', and the 
above mentioned Vp^s is its distribution. 

We use the notation 0{p, 8) for the probability that, in the final config- 
uration, O is in an infinite occupied cluster: 

e{p,5) :=Pp,a(0^oo). 

Note that O is occupied in the final configuration if and only if the above 
mentioned enhancement was successful, or O belonged initially (before the 
catastrophe) to a non-empty but finite occupied cluster. This gives 

Vp,s{0 is occupied ) = 5 + {I - 6){p - 9{p)). 
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Also note that, in the case that p < Pc, nothing happens in the above 
catastrophe, so that in the final configuration the sites are independently 
occupied with probability p+ {1 — p)5. Formally, li p < Pc, then 



where we use the notation Vp for the product measure with parameter p. In 
particular. 



for each S > 0. 

Remark 1.1. Most of what we said above has straightforward analogs for 
arbitrary countable graphs, hut there are subtle differences. For instance, on 
the cubic lattice it has not yet been proved that 9{pc) = (although this is 
generally believed to be true). So, for that lattice, (2) with p = Pc, and hence 
(3), are not rigorously known. 

It is also clear from the construction that Pp^s stochastically dominates Ps- 
Hence, if S > Pc then 9{p, 6) > 9{S) > for all p. 

It turns out (see Proposition 3.1 of [1]) that, ii p > Pc, a 'non-negligible' en- 
hancement is needed after the catastrophe to create again an infinite occu- 
pied cluster. More precisely, for each p > Pc there is a 5 > with 6{p, 5) = 0. 
A much more difficult question is whether the needed enhancement goes to 
as p J, Pc- By (3) one might be tempted to reason intuitively that this is 
indeed the case. In [1] it was shown that for the analogous model on the 
binary tree this is correct. However, in [1] a conjecture is presented which 
says, in particular, that for the square lattice (and other planar lattices) 
there is a (5 > for which 6{p, S) = for all p > Pc- In Section 4 of [1] and 
in [2] we showed remarkable consequences for certain forest-fire models. 

Note that, since 6{pc, 5) > 0, the above conjecture says that the function 
9{-, •) has discontinuities at points of the form {pc, 5) with 6 sufficiently small. 
This naturally raises the question whether this function is continuous in the 
complement of a region of such form: is there a S > such that 9{-, •) is 
continuous outside the set {pc} x [0, 6]? In the next subsection we state that 
this is indeed the case, and give a summary of the methods and intermediate 
results used in the proof. At the end of Section 6 we point out why our proof 
does not work at points {pc, 5) with small S. We hope our arguments provide 
a better understanding of the earlier mentioned conjecture and will trigger 
new attempts to prove (or disprove) it. 



'PptS — 'Pp+{l—p)5 1 



(2) 



9{pc,6)=9{pc + {l-Pc)S)>0, 



(3) 
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1.2 Outline of results 



The conjecture mentioned in the previous subsection raises the natural ques- 
tion whether ^(•, •) is continuous outside the indicated 'suspected' region. 
The following theorem states that this is indeed the case. 

Theorem 1.2. There is a 5 & (0, 1) such that the function 9{-, •) is contin- 
uous outside the segment {pc} x {0,6). 

As could be expected, the proof widely uses tools and results from 
ordinary percolation. However, the dependencies introduced by the self- 
destructive mechanism cause complications. Until recently, a serious obsta- 
cle was the absence of a suitable RSW-like theorem. This obstacle could be 
removed by the use of (a modified and somewhat stronger form of) a recent 
theorem of BoUobas and Riordan ([4]). 

A rough outline of the proof of Theorem 1.2, and the needed intermediate 
results that are interesting in themselves, is as follows: In section 2 we list 
some basic properties of our model, which will be used later. The results 
in Section 3, which are also contained in the recent PhD thesis [6] of one 
of us, show that if 6{-, •) is strictly positive in some open region, then it is 
continuous on this region. It is also shown that if 9{p,S) = 0, then 6{-, •) 
is continuous at {p,d). These two results reduce the proof of Theorem 1.2 
to showing that if 9{p,6) > and p ^ pc, then 9{p,5) > in an open 
neighborhood of {p,S). This in turn requires a suitable finite-size criterion 
(see below) for sdp. In Section 4 we give the modified form of the BoUobas- 
Riordan theorem. This is used in Section 5 to obtain the above mentioned 
finite-size criterion. Finally, in Section 6 we combine these results and prove 
the main theorem. 

We end this section with the following remark: When we say that a 
function / is 'increasing' ('decreasing') this should, unless this is preceded 
by the word 'strictly', be interpreted in the weak sense: x < y implies 
fix) < f{y). 

2 Basic properties 

In this section we state some basic properties which will be used later. 

First some more terminology and notation: li v = (^1,^2) and w = 
{wi,W2) are two vertices, we let \v — w\ denote their (graph) distance — 
wi \ + \v2 — W2\- By B{v, k) and dB{v, k) we denote the set of vertices w for 
which — is at most k, respectively equal to k. For V,W C 1? , we define 
the distance between V and W as min{|i; — w\ : v &V,w € W}. 
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Recall that Vp^s denotes the sdp distribution (that is, the distribution of 
the collection (Zj, i G W?) defined in Subsection 1.1). This is a distribution 
on 17 := {0, 1}^ (with the usual cr- field). Elements of are typically denoted 
by uj{= {iOi, i G Z^)), cr etc. We write w < a if < (Tj for all i. 
Let F be a set of vertices and A an event. We say that A lives on V if 
u ^ A and cti = uoi for all z G V , implies a ^ A. And we say that A is a 
cylinder event if A lives on some finite set of vertices. As usual, wc say that 
A is increasing if a; G ^4 and < ai for all i, implies a £ A. The first two 
lemma's below come from Section 2.2 and 2.4 respectively in [1]. 

Lemma 2.1. Let A and B he two increasing cylinder events. We have 
rp,s{AnB)>Pp,s{A)Vp,5{B). 



As to monotonicity, it is obvious that the sdp model has monotonicity 
in 6: If Si > 62, then Pp,Si stochastically dominates Vp,S2- Although there 
seems to be no 'nice' monotonicity in p we have the following result. 

Lemma 2.2. Ifp2>Pi and p2 + {1 — P2)S2 < Pi + — Pi)Si, then 

'Ppi,5i dominates Vp^^s^- 



The next result is about 'almost independence' of cylinder events which 
live on widely separated sets. As usual, the lattice which has the same ver- 
tices as the square lattice but where each vertex has, besides the four edges 
to its nearest neighbours, also four 'diagonal edges' is called the matching 
lattice (of the square lattice) . To distinguish paths and circuits in the match- 
ing lattice from those in the square lattice, we use the terminology *-paths 
and *-circuits. 

Lemma 2.3. Let k be a positive integer and let V and W be subsets of 
that have distance larger than 2k. Further, let A and B be events which live 
on V and W respectively. Then 

\Vp4AnB)-Vp4A)Vp4B)\ < (4) 
2{\V\ + \W\)Vp (3 vacant *-circuit surrounding O and some vertex in dB{0,k)) . 
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Proof. Recall how we formally defined the sdp model in terms of random 

variables X, Y and Z. Wc use a modification of those variables: Let X and 
Y be as before, but in addition to X* and Z we now define and Z^f^ 

by 

*(/) _ J 1 if = 1 and there is no X- occupied path from i to dB{i, k) 
otherwise; 



yif) _ v*{f) 



X^'^'yYi. (5) 

Let Vp'^^ denote the distribution of Z^-f ^ . It is clear that the random variables 

(f) ' if] 
Zl ,i &V are independent of the random variables Z- ,i gW, and hence 

vi/JiAnB) = vifJ{A)vl^J{B). (6) 

( f) 

Also note that if Zi ^ -^j , then the X— occupied cluster of i intersects 
k) but is finite. Hence there is an X— vacant circuit in the matching 
lattice that surrounds i and some site in dB{i,k). Hence, since the X— 
variables are Bernoulli random variables with parameter p, we have for any 
finite set K of vertices and any event E living on K, 



- VI{1{E)\ < P{Zk + Zi^) < (7) 
\K\Vp (3 a vacant *-circuit surrounding O and some vertex in dB{0, k)) . 

The lemma now follows easily from (6) and (7) □ 

Our last result in this section is on the uniqueness of the infinite cluster. 

Lemma 2.4. If9{p,S) > 0, then 

Vp^si^ O' unique infinite occupied cluster ) = 1. 



Proof. Prom the earlier construction of the sdp model in terms of the X— 
and Y variables, it is clear that Vp^s is stationary and ergodic. It is also clear 
that in the sdp model the conditional probability that a given site is occupied 
given the configuration at all other sites, is at least 5. So this model has 
the so-called positive finite energy property. The result now follows from an 
extension in [9] of the well-known Burton-Keane ([7]) uniqueness result. □ 
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3 Partial continuity results 



In this section we first prove that in the sdp model the probabihties of 
cylinder events are continuous functions of (p, 5) . Next we prove that the 
function ^(•,-) is continuous at (p, (5) if 9{p,6) = or there is an open 

neighborhood of (p, S) on which 9 is strictly positive. Note that, once we have 
this, the proof of Theorem 1.2 is basically reduced to showing that \ip ^ Pc 
and 9{pc,S) > 0, then 9{-, •) is strictly positive on an open neighborhood of 
{P,S). 

Lemma 3.1. Let A be a cylinder event. The function {p,S) — >■ Vp^siA) is 
continuous on [0,1]^. 

Remark 3.2. The proof (see below) uses the well-known fact that 9{pc) = 0. 
For many lattices (e.g. the cubic lattice) this fact has not been proved. For 
those lattices the arguments below show that the function in the statement 
of 3.1 is continuous on [0, 1]^ \ {{pc} x [0, 1]). 

Proof. Let A be an event which lives on some finite set V. Recall the con- 
struction of the sdp model in terms of random variables X, Y and Z. Let, 
for a E ft, ay denote the tuple {ai,i G V). It is clear that the distribution 
of Xy is a function of p only, and that, conditioned on Xy, the probability 
that Zy (z A is a polynomial (of degree \V\) in 5. Therefore it is sufficient 
to prove that, for each a G {0,1}^, the function f '■ p —>■ V{Xy = a) 
is continuous. Recall that the X— variables are Bernoulli random variables 
(with parameter p). Now let < pi < P2- In a standard way, by introducing 
independent, uniformly on the interval (0, 1) distributed random variables 
Ui,i & I? , Vie can suitably couple two collections of Bernoulli random vari- 
ables with parameters pi , respectively p2 . Such argument easily gives that 
\fip2) — f{pi)\ is less than or equal to the sum over z G F of 

'P{Ui € (pi,P2)) +V{i is in an infinite p2-open but not in an infinite pi-open cluster), 
which equals 

\v\{p2-Pi + e{p2)-e{pi)). 

The lemma now follows from the continuity of ^(.). □ 

Proposition 3.3. Let {p,S) G [0,1]^. If (a) or (b) below holds, the function 

9{-,-) is continuous at {p,S). 

(a) 9{-, ■) > on an open neighborhood of ip,6). 

(b) 9{p,5) = 0, 
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Proof. For this (and some other) results it is convenient to describe the sdp 
model in terms of Poisson processes: Assign to each site, independently 
of the other sites, a Poisson clock with rate 1. These clocks govern the 
following time evolution: Initially each site is vacant. Whenever the clock 
of a site rings, the site becomes occupied. (If it was already occupied, the 
ring is ignored). Note that if occupied sites would always remain occupied, 
then for each time t, the configuration at time t would be a collection of 
independent Bernoulli random variables with parameter 1 — exp(— i). In 
particular, before and at time tc, defined by the relation Pc = 1 — exp(— tc), 
there would be no infinite occupied cluster, but after tc there would be 
a (unique) infinite cluster. However, we do allow occupied sites to become 
vacant, although only once, as follows: Fix a time r, a parameter of the time 
evolution. At time r all sites in the infinite occupied cluster become vacant. 
(If there is no infinite occupied cluster, which is a.s. the case if r < tc, 
nothing happens). After time r we let the evolution behave as before; that 
is, each vacant site becomes occupied when its Poisson clock rings. Let, for 
this time evolution with parameter r, VT,t denote the distribution of the 
configuration at time t , and let 

9{T,t) = Vr,t{0 is in an infinite occupied cluster ). (8) 

It is easy to see that 

Vr,t = Vp,S, (9) 

where p = 1 — exp(— r) and 6 = 1 — exp(— (t — r)). It is also easy to see 
that Pr,t is stochastically decreasing in r and stochastically increasing in t. 
In fact this is the key behind Lemma 2.2. 

Now we come back to the proof of Proposition 3.3. From (8) and (9) 
we get (since the map between pairs {p, 6) and (r, t) in (9) is continuous) 
that this proposition is equivalent to saying that if 9{T,t) = or t tc and 
6 is strictly positive on an open neighborhood of {T,tc), then 9 is continu- 
ous at (r, t). To prove this equivalent form of Proposition 3.3 we use ideas 
from [3] . The introduction of pairs (r, t) as replacement of (p, 6) not only 
has the advantage that, as we already saw, we now have a more suitable 
form of monotonicity, but, more importantly, that we now have a more 'de- 
tailed' structure (the Poisson processes) in the background which gives the 
appropriate 'room' needed to get a suitable modification of the arguments 
in [3]. 

Let (r, t) be as above. Divide the parameter space in four 'quadrants', 
numbered / to IV: 
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/:= [0,r] X [t,oo), 

II := [r, oo) X [t, oo), 

III := [r,oo) X [0,t], 

IV := [0,t] X [O.t]. 

Note that it is sufficient to prove that for each monotone sequence (tj, ii)i>o 
that hes in one of the above quadrants and converges to (r, t) , one has 

hm 9{Ti,ti) = 9{T,t). 

i—^oo 

We handle each of the quadrants separately. 

Quadrant I) This is easy and corresponds to the (easy) proof of right con- 
tinuity of ordinary percolation: Let (rj) be a monotone sequence which 
converges from below to r and let (tj) be a monotone sequence which con- 
verges from above to t. Let An denote the event that there is an occupied 
path from O to dB{0,n). By monotonicity and Lemma 3.1 we have that 

For each i,Vr,,ti{An) j Oin^ti) as n ^ oo; (10) 
VrMn) i 0{t, t)asn^ oo; (11) 
For each n,VTi,ti{An) i Vr,t{An) as i ^ oo, (12) 

From these three statements it is easy to see that 6{Ti, ti) tends to 6{t, t) as 
z ^ oo. 

Quadrant III) Let (r^) be a monotone sequence which converges from 
above to r and {ti) a monotone sequence which converges from below to t. 
By the earlier monotonicity arguments, the sequence 6{Ti,ti) is increasing 
in i, and has a limit smaller than or equal to 0{T,t). So for the situation 
where 9{T,t) = 0, the proof is done. Now we handle the other situation: we 
assume 6 is positive in an open neighborhood of (t, t). For this situation 
considerable work has to be done. Note that in the dynamic description 
given earlier in this section, the underlying Poisson processes were the same 
for each choice of the model parameter r. This allows us (and we already 
used this to derive some monotonicity properties) to couple the models with 
the different r^'s and r. 

Let, for s < u, Cg^u denote the occupied cluster of site O at time u 
in the process with parameter s (that is, under the time evolution where 
the infinite occupied cluster is destroyed at time s). Further, we use the 
notation lo{s,u) for the configuration at time u in that model. It is also 



9 



convenient to consider uj{u), the configuration at time u in the model where 
no destruction takes place. (So, lo„{u),v £ 7?, are independent — 1 valued 
random variables, each being 1 with probability 1 — exp(— n)). Again we 
emphasize that all these models are defined in terms of the same Poisson 
processes. Prom monotonicity (note that C^-.^t. C C^-.^^^t.^^ for all i) it is 
clear that 

lim e{n,ti) = p{3i\Cr,,t,\ = ^), 

and 

9{T,t) - lim e{n,U) = V{\Cr,t\ = oo, Vi \Cr,,u\ < oo). (13) 
i—*oo 

So we have to show that the r.h.s. of (13) is 0. Pix a j with the property 

that 9{Tj,tj) > 0. Such j exists by the condition we assumed for (r, f). To 
show that the r.h.s. of (13) is 0, it is sufficient (and necessary) to prove the 
following claim: 

Claim 

Apart from an event of probability 0, the event {\Cr,t\ = oo} is contained in 
the event that there is a k > j for which |Crfc,tfcl = 

So suppose \Cr,t\ = oo. By our choice of j we may assume that Lu{Tj,tj) has 
an infinite occupied cluster, and by Lemma 2.4 that this cluster is unique. 
We denote it by Ij. If O G Ij we are done. Prom monotonicity and the 
uniqueness of the infinite cluster (see Lemma 2.4), we have Ij C Cr^t- Hence 
there is a finite path vr from O to some site in Ij such that u}{t, t) = 1 on vr. 
Since, a.s. there are no vertices whose clock rings exactly at time t or r, we 
may assume that for every site v on tt, (a) or (b) below holds: 

(a) The clock of v rings in the interval (r, t) . 

(b) u;v{t) = 1 but the occupied cluster of v in u!{t) is finite. 

If (a) occurs we define: 

:= min{i : i > j and the clock of v rings in (r^, tj)}. 

Note that then, by the monotonicity of the sequence (Ti,ti), the clock of v 
rings in the interval (r;,t;) for all I > iv If (a) does not occur, (b) occurs, 
and hence there is a finite set of sites on which uj{t) = and which 
separates v from oo. Then we use the following alternative definition of i^: 

iv ■= min{z : i > j and a;(rj) = on K^}. 

This minimum exists since is finite and (again) we assume that no Poisson 
clock rings exactly at time r. Now let 
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k := maxiy, 

which exists since tt is finite. 

From the above procedure it is clear that uj^{Tk,tk) = 1 for all v on tt. 
Further, since k > j and by monotonicity, also ujy{Tk,tk) = 1 for all v Ij. 
Since tt is a path from O to Ij this implies that Ij is contained in C^-^ and 
hence that |CT-j.,tj,| = oo. This proves the Claim above. 

Quadrants II) and IV) 

The required results for these quadrants follow very easily from monotonicity 
and the above results for quadrants I and III: Let (Tj,tj) be a sequence in 
quadrant II that converges to (r, t). We have, by earlier stated monotonicity 
properties, 

o{ri,t)<e{n,ti)<Hr,ti)- 

Since the sequence {Ti,t)) lies in quadrant III and the sequence {T,ti) 
lies in quadrant I, the upper and lower bound both converge to 6{t, t). This 
completes the treatment of quadrant II. Quadrant IV is treated in the same 
way. This completes the proof of Proposition 3.3 □ 

4 An RSW-type result 

For our main result we need to prove that if the crossing probability of an 
n by n square goes to 1 as n — > oo, then also the crossing probability of 
a (say) 3n x n rectangle in the 'difficult direction' goes to 1 as n ^ oo. 
Such (and stronger) results were proved for ordinary percolation in the late 
nineteen seventies by Russo, and by Seymour and Welsh, and therefore 
became known as RSW theorems. Their proofs used careful conditioning 
on the lowest horizontal crossing in a rectangle, after which the area above 
that crossing was treated, and a new, vertical crossing in that area was 
'constructed'. Such arguments work for ordinary percolation because there 
the above mentioned area can be treated as 'fresh' territory. However, they 
usually break down in situations where we have dependencies, as in the sdp 
model. 

Recently, BoUobas and Riordan ([4]) made significant progress on these 
matters. For the so-called Voronoi percolation model they proved an RSW 
type result. That result is one of the main ingredients in their proof that 
the critical probability for Voronoi percolation equals 1/2 (which had been 
conjectured but stayed open for a long time). Although they explicitly 
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proved their RSW type result only for the Voronoi model, their proof works 

(as they remark in their paper) for a large class of models. The result we 
needed is a little stronger than that of [4] . The rest of this section is organised 
as follows. First we give a short introduction to Voronoi percolation. Then 
we state the above mentioned RSW-like theorem of [4] , and point out where 
and how its proof needs to be modified to obtain the stronger version. Finally 
we state the analog for the sdp model and explain why the proof for the 
Voronoi model works for this model as well. 

4.1 The Voronoi percolation model 

We start with a brief description of the Voronoi percolation model. The 
(random) Voronoi percolation model is as follows: Let Z denote the (ran- 
dom) set of points in a Poisson point process with density 1 in the plane. 
This set gives rise to a random Voronoi tessellation of the plane: Assign 
to each z ^ Z the set of all a; G for which z is the nearest point in Z. 

The closure of this set is called the (Voronoi) cell oi z. It is known that 
(with probability 1) each Voronoi cell is a convex polygon, and that two 
cells are either disjoint or share an entire edge. In the latter case the two 
cells are said to be neighbours or adjacent. This notion of adjacency gives, 
in a natural way, rise to the notion of paths, clusters etc. 

Now consider the percolation model where each cell, independently of 
everything else, is coloured black with probability p and white with prob- 
ability \ — J). Based on analogies with ordinary percolation (in particular 
with the self-matching property of the usual triangular lattice) it has been 
conjectured for a long time that the critical value for this percolation model 
is 1/2: for p < 1/2 there is (a.s.) no infinite black cluster, but for p > 1/2 
there is an infinite black cluster (a.s.). As we said before, this was recently 
proved rigorously by Bollobas and Riordan ([4]), and a key ingredient in 
their proof is an ingenious RSW-like result. 

4.2 The RSW-like result for Voronoi percolation 

As in [4] we define, for the Voronoi percolation model with parameter p, 
/p(p, s) as the probability that there is a horizontal black crossing of the 
rectangle [0, ps\ x [0, s] . The following is Theorem 12 in [4] 

Theorem 4.1. (Bollobas and Riordan) Let < p < 1 be fixed. 
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If limmf/p(l,s) >0, (14) 

S—i-OO 

then limsup /p(p, s) > for all p > 0. 

s—*oo 

Studying the proof we realised that the condition can be weakened, so 
that the following theorem is obtained: 

Theorem 4.2. Let < p < 1 be fixed. 

If lim sup fp{p, s) > for some p > 0, 

s-^oo 

then lim sup fp{p, s) > for all p > 0. 

s— »oo 



As we shall point out, this somewhat stronger Theorem 4.2 can be proved in 
almost the same way as Theorem 4.1. But see Remark 4.4 about the global 
structure of the proof. First note that Theorem 4.2 is (trivially) equivalent 
to the following: 

Theorem 4.3. Let < p < 1 be fixed. 

If limsup fp{p, s) = for some p > 0, (15) 

s— ►oo 

then 

limsup/p(p, s) = for all p > 0. (16) 

This is the form we will prove, following (with some small changes) the 
steps in [4]. 

Proof. (Theorem 4.2 and 4.3). Since p is fixed we will omit it from our 
notation. In particular wc will write / instead of fp. 

First we will rewrite the condition (15) in Theorem 4.3: If limsupg^g^ f{p, s) = 
for some p < 1 then, since f{p,s) is decreasing in p, this limsup is for 
all p > 1. Moreover, the well-known pasting techniques from ordinary per- 
colation show easily that if limsup^^oo f{p, s) > for some p > 1, then this 
limsup is positive for all p' > p, and hence (using again monotonicity of / 
in p) for all p > 1. Equivalently, if limsupg^^ f {p, s) = for some p > 1, 
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then this hmit equals for all p > 1. Hence, the condition in Theorem 4.3 
is equivalent to 

limsup/p(p, s) = for all p > I, (17) 
and this is also equivalent to condition (3) in Section 4 of [4]: 

limsup/p(/3, s) = for some p > 1, (18) 

s— >oo 

Wc will assume (17) (or its equivalent form (18)) and show how, follow- 
ing basically the proof of Theorem 4.1, the equation in (16) can be derived 
from it for all p > 1/2. Then we make clear that, for each k, very similar ar- 
guments work for 1/k instead of 1/2, which completes the proof of Theorem 
4.3. 

Remark 4.4. BoUobds and Riordan prove their theorem by contradiction: 
They assume (as we do here) (18) above, and, m,oreover they assume (14) 
(equation (2) in Section 4 of their paper). Then, after a number of steps 
(claims), they reach a contradiction, which completes the proof. However, 
most of these steps do not use the 'additional' assumption (14) at all. We 
found a 'direct' (that is, not by contradiction) proof, as sketched below, more 
clarifying since it leads more easily to further improvements. For our goal 
most of the steps (claims) in the proof in [4] remain practically unchanged. 
Therefore we (re) write only some of them in more detail (Claim 1 is stated 
to give an impression of the start of the proof, and Claim 4 because that 
already gives a good indication of the strong consequences of (18) ). For 
the other claims we only describe which changes have to he made for our 
purpose. 

First some notation and terminology: Tg is defined as the strip [0, s] X 
M. An event is said to hold with high probability, abbreviated whp if its 
probability goes to 1 as s — > cx3 (and all other parameters, e.g. p and e are 
fixed). 

Claim 1 (Claim 12.1 in [4]). 

Let e > 6e fixed, and let L be the line-segment {0} x [— es,es]. Assuming 
that (18) holds, the probability that there is a black path P in Tg starting 
from L and going outside S' = [0, s] x [—(1/2 + 2e)s, (1/2 + 2£)s] tends to 
as s ^ DO. 
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This claim is exactly the same as in [4] , except that in their formulation not 
only (18) but also (14) is assumed. However, their proof of this claim does 
not use the latter assumption. 

The above, quite innocent looking claim, leads step by step to stronger 
and eventually very strong claims. We will not rewrite Claim 2 and Claim 
3; like Claim 1, they are exactly the same as their corresponding Claims 
(12.2 and 12.3 respectively) in [4], except that we do not assume (14). And, 
again, the proof remains as in [4]. 

Claim 4 (Claim 12.4 in [4]) Let C > be fixed, and let R = Rg be the s 
by 2Cs rectangle [0, s] x [-Cs,Cs]. For < j < 4, set Rj = [js/100, (j + 
96)s/100] X [—Cs,Cs]. Assuming that (18) holds, whp every black path P 
crossing R horizontally contains 16 disjoint black paths Pi, 1 < i < 16, 
where each Pi crosses some Rj horizontally. 

Again, in the formulation in [4] also (14) is assumed, but this is not used in 
the proof. Following [4] we now define, for a rectangle R, the random variable 
L{R) as the minimum length of a black path crossing R horizontally. (More 
precisely, it is the minimum length of a piecewise-linear black curve that 
crosses R horizontally). If there is no horizontal black crossing of R we 
take L{R) = oo. A complicating property of L is that if Ri and R2 are 
two disjoint rectangles, L{Ri) and L(i?2) are not independent (no matter 
how large the distance between the two rectangles). Therefore, below Claim 
12.4 in their paper, Bollobas and Riordan introduce a suitable modification 
L. The key idea is that whp the colours inside a rectangle with length and 
width of order s, are completely determined by the Poisson points within 
distance of order o(s) of the rectangle. 

There are many suitable choices of L, and we will not rewrite the precise 
definition given in [4], but only highlight the following three key properties 
(which neither use (18) nor (14)): 

L{Rs) = L{Rs), whp, (19) 
L{Rs) > s, whp, (20) 

and: 

Claim 5. (Claim 12.5 in [4]). Let R\ and R2 be two s by 2s rectangles, 
separated by a distance of at least s/100. // s is large enough, then the 
random variables L{Ri) and L{R2) are independent. 

Remark 4.5. In fact, the independence property of L is only used in the 
proof of Claim 12.6 in [4], and there it could be replaced by the following 
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property (which follows from (19) and Claim 5): 

For each e > there is au such that for all s > u and all s by 2s rectangles 
Ri and R2 that are separated by a distance of at least s/100, we have 

sup \PiL{Ri) < X, L{R2) <y) - P{L{Ri) < x)P{L{R2) < y)\ < e. 

x,y>0 

Now choose an arbitrary number fj smaller than 10~^. This deviates from 
the choice of rj by Bollobas and Riordan, who add an extra condition, related 
to their assumption of (14). Define 

i{s) = sup{s : V{L{Rs) < x) < f]}. (21) 

This definition of t is the same in form as that of t in [4] (see two 
lines below (16) in [4]); however our way of choosing fj was different. A 
consequence of this difference is that, in our setup, t{s) can be 00. As in [4], 
we do have that 

i{s) > s, for all s. (22) 

Claim 6. (Claim 12.6 in [4]). 
Let Rs be a fixed 0.96 s by 2s rectangle. If (18) holds, then 

V [LiRs) < i{0.A7s)) < 200 f. 

This statement is the same as in [4], except that in [4] also (14) is as- 
sumed, and that we use t and fj instead of t, respectively rj. The proof is 
the same as in [4]. 

Prom the above (in particular Claim 4, Claim 6 and (19)), the following 
quite startling Proposition (which, essentially is equation (18) in [4]) now 
follows quite easily. 

Proposition 4.6. (Corresponds with (18) in [4]). If (18) holds, then, for 
all sufficiently large s, 

i{s) > 16t(0.47s). (23) 
Proof. Practically the same as the proof of equation (18) in [4]. □ 
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Now Theorem 4.3 follows in a few lines from this proposition and the 
definition of t{s): It is easy to show (see the arguments below (18) in [4]) that 
(22) and Proposition 4.6 together imply that t{s) > s^, for all sufficiently 
large s. Hence, by the definition of t{s) (and by (19)) we get that, for all 
sufficiently large s, 

V{L{Rs) < s^) < r). (24) 

It is also easy to show (see the arguments below equation (19) in [4]) 
that 

7^(5^ < L{Rs) < oo) ^ as s ^ oo. (25) 

Combining (24) and (25), and recalling that L{Rs) = oo iff there is no 
horizontal black crossing of Rg, immediately gives 

limsupP(3 a horizontal black crossing oi Rg) < fj. 

S—KX> 

Now, since r) was an arbitrary number between and 10~^, we get 
lims_^ooP(3 horizontal black crossing of Rg) = 0, that is, 

/(l/2,s) ^0, as s ^ cx). (26) 

Note that in the last part of the above arguments (after Claim 4) we 
worked in particular with s by 2s rectangles. A careful look at the arguments 
shows that the choice of this factor 2 is, in fact, immaterial: if we would 
take s by 3s rectangles or, more generally, fix an > 2 and take s by Ns 
rectangles, the arguments remain practically the same. To see this, one can 
easily check that in Claims 1-4 (Claims 12.1 - 12.4 in [4]) the factor 2 plays 
no role at all: here the rectangles under consideration are s by 2Cs, where 
C is a fixed but arbitrary positive number. Further, the proof of Claim 5 
remains the same when, for some fixed positive number C, we replace the 
factor 2 by 2C. And, the definition of i{s) (see (21)), which was given in 
terms of s by 2s rectangles, has, for each C > 0, an obvious analog for s by 
2Cs rectangles: 

tcis) := sup{s : V{L{R^) < x) < r}}, (27) 

where, for each s, i?^ is some fixed s by 2Cs reactangle. 

In the generalization of Claim 6 (Claim 12.6 in [4]) we now fix C > 1, 
and take R^ := [0, 0.96s] x [—Cs, Cs]. In the proof of this Claim we have to 
replace, on the vertical scale, s by Cs. For instance, the segments Lj, which 
in the original proof in [4] have length 0.02s, will now have length 0.02Cs, 
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and Rq and Ri which in the original proof are 0.47s by 2 x 0.47s rectangles, 
are now 0.47s by 2C0.47s rectangles. In this way we get if (18) holds, for 
each fixed C > 1 the following analog of (26): 

/(1/(2C), s) ^0, as s ^ oo. (28) 
This proves Theorem 4.3 and hence Theorem 4.2. 

□ 

In the above we were dealing with black horizontal crossings. Obviously, 
completely analogoTis results hold for white horizontal crossings: If we de- 
note (for a fixed value of the parameter p of the Voronoi percolation model), 
the probability of a vertical white crossing of a given ps by s rectangle by 
g{p, s), we have that if lim^-joo dip, s) = for some p > 0, then this limit is 
for all p > 0. Since a rectangle has either a horizontal black crossing or a 
vertical white crossing (and hence g{p, s) = 1 — f{p, s)) this gives: 

Corollary 4.7. 

// lim f{p, s) = 1 for some p > 0, (29) 

s— »oo 

then lim f{p, s) = 1 for all p > 0, 



4.3 An RSW analog for self-destructive percolation 

In the previous subsection we considered (and somewhat strengthened) an 
RSW- like result of Bollobas and Riordan ([4]) for the Voronoi percolation 
model. Only a few properties of the model are used in its proof. As remarked 
in [4] (at the end of Section 4; see also [5], Section 5.1), these properties 
are basically the following: First of all, crossings of rectangles are defined 
in terms of 'geometric paths' in such a way that, for example, horizontal 
and vertical black crossings meet, which enables to form longer paths by 
pasting together several small paths. Further, a form of FKG is used (e.g. 
that events of the form 'there is a black path from A to are positively 
correlated. Also some symmetry is needed. Bollobas and Riordan say that 
"invariance of the model under the symmetries of suffices, as we need 
only consider rectangles with integer coordinates". Finally, some form of 
asymptotic independence is needed (see Remark 4.5). Similar considerations 
hold wrt the somewhat stronger Theorem 4.2. 

Using the results in Section 2, is not difficult to see that the sdp model has 
the above mentioned properties: 
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• The indicated geometric properties are just the well-known intersection 
properties of paths in the square lattice (and in its matching lattice). 

• Lemma 2.1 gives the needed FKG-like properties. 

• Asymptotic independence: Note that for p < Pc the sdp model is an 
ordinary percolation model, where this property is trivial. If p > Pc, 
then 1 — p is smaller than the critical probability of the matching 
lattice. In that case the needed asymptotic independence (of the form 
described in Lemma 4.5) comes from Lemma 2.3 and the well-known 
exponential decay theorems for ordinary subcritical percolation. 

• The sdp model on the square lattice clearly has all the symmetries of 

Further, to carry out for the sdp model the analog of the arguments that 
led from Theorem 4.3 to Corollary 4.7, we note that the random collection 
of vacant sites on the matching lattice clearly also has the above mentioned 
properties. So we get the following theorem for the sdp model: 

Theorem 4.8. The analogs of Theorems and 4-3 and Corollary 4 -"^ hold 
for the self-destructive percolation model. In particular, let for the sdp model 
with parameters p and 6, f{p, s) = fp,s{p, s) denote the probability that there 
is an occupied horizontal crossing of a given ps x s rectangle. We have 

If lim f(p, s) = 1 for some p > 0, (30) 

s— »oo 

then lim f{p, s) = 1 for all p > 0. 

s— >oo 

In the next sections this result will play an important role in the com- 
pletion of the proof of Theorem 1.2. In particular, in Section 5 it will be 
used to prove a finite-size criterion for supercriticality of the sdp model. 

5 A finite-size criterion 

The main result of this section is a suitable finite-size criterion for super- 
criticality of the sdp model. The overall structure of the argument is similar 
to that in ordinary percolation (see [8]), but the dependencies in the model 
require extra attention. One of the main ingredients, a suitable RSW-like 
theorem for this model, was obtained in the previous section. 
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Theorem 5.1. Let f = fp^s in Theorem 4.8. There is a universal con- 
stant a > and there is a decreasing function N : (pc, 1) ^ N such that for 
all p > Pc and all 6 > the following two assertions, (i) and (ii) below, are 
equivalent. 

i. e{p,6) > 0. (31) 
ii3n > N{p) such that /p,j(3, n) > 1 — a. 

Remark 5.2. In ordinary percolation N{p) can be taken constant 1. Remark 
6. 1 below explains the impact of this difference. 

Proof. Consider for each n € N the events 

^4 = {3 a vertical vacant *-crossing of [0, 9n] x [0,3n]}; (32) 
S = {3 a vertical vacant *-crossing of [0, 9n] x [0,n]}; 
C = {3 a vertical vacant *-crossing of [0, 9n] x [2n, 3n]}. 

Let h{p, n) denote the probability of a vertical vacant crossing (in the 
matching lattice) of a pn by n box. So, h{p,n) = 1 — f{p,n). Clearly, 
Vp^s{B) = Vp^s{C) = ^p,5(9,ra) and Vp^5{A) = /ip5(3, 3n). It is also clear 
that A C B D C. Prom this. Lemma 2.3, the fact that the r.h.s. of (4) 
is decreasing, and the well-known exponential decay results for ordinary 
subcritical percolation applied to (4), it follows that there is an increasing, 
function cf) : {pc, 1) —>■ {0, 00) such that for all p > Pc, 

h{3, 3n) < h{9, nf + exp(-n^(p)). (33) 

Further note that if the event B occurs, there must be a vacant vertical *- 
crossing of one of the rectangles [0, 2>n] x [0, n] , [2n, 5n] x [0, n] , [4n, 7n] x [0, n] , 
[6n, 9n] x [0, n], or a vacant horizontal *-crossing of one of the rectangles 
[2n,3n] x [0,n], [4n, 5n] x [0,n], [6n, 7n] x [0,n]. 

Hence 

/i(9, n) < 4/1(3, n) + 3/1(1, n) < 7/i(3, n), (34) 
which combined with (33) gives 

/i(3, 3n) < 49/i(3, nf + exp(-n<^(p)). (35) 
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Take a so small that 49a^ < q;/4. Let, for each ■p > Pc, N{p) be the 
smallest positive integer for which 

exp(-7V(p)(?!)(p)) < a/4. 

Sine (j) is increasing, N is decreasing in p. 

Now suppose p > Pc and 6 G (0, 1) are given and suppose that (ii) holds. So 
there exists an n that satisfies: 

exp(— n(^(p)) < a/4 and h{3,n) < a. (36) 

Prom (35) , (36) and the choice of a we get 

h{3, 3n) < 49a2 + a/4 < a/4 + a/4 = a/2, (37) 

and 

exp(-3n(/.(p)) < {a/4f < (a/2)/4. 

Hence, (36) with n replaced by 3n, and a replaced by a/2 holds. So we 
can iterate (37) and conclude that, for all integers A: > 0, 
h{3, 3*^n) < a/(2'=). 

The last part of the argument is exactly as for ordinary percolation: 
Note that if none of the reactangles [0,3^*^+^7^] x and [0, S^'^+^n] x 

[0, 32'=+2n], fc = 0, 1, 2, • • • has a white * -crossing in the 'easy' (short) direc- 
tion, then each of these rectangles has a black crossing in the long direction. 
Moreover, all these black crossings together form an infinite occupied path. 
Hence, 

oo oo 

9{p,6) > 1 -^/i(3,n3*^) > 1 -a^(l/2)'' = 1 - 2a > 0. 

k=0 k=0 

This proves that (ii) implies (i). 

Now we show that (i) implies (ii): Suppose 0{p,S) > 0. Then there is (a.s.) 
an infinite occupied cluster, and by Lemma 2.4 this cluster is unique. From 
the usual spatial symmetries, positive association, and the above mentioned 
uniqueness one can, in exactly the same way as for ordinary percolation (see 
[10], Theorem 8.97) show that /(l,n) ^ 1 as n ^ oo. By Theorem 4.8 it 
follows that also /(3, n) ^ 1 as n — > oo; so (ii) holds. □ 
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6 Proof of Theorem 1.2 



We are now ready to prove Theorem 1.2: 

Proof. For p < pc, we have (see Section 1) 9{p,6) = 9{p + (1 — p)S), so 
that continuity fohows from continuity for ordinary percolation, li p = pc 
and 6 > Pc + s foT some £ > 0, then (triviahy) there is a neighborhood of 
{p, 6) where the sdp model dominates ordinary percolation with parameter 
Pc + e/2 > Pc] hence 6{-, •) > on this neighborhood, and Proposition 3.3 
implies continuity of 9{-, •) at {pc, 5). (In fact, by combining this argument 
with an Aizenman-Grimmett type argument, one can extend this result 
and show that there is an e > such that 9{-,-) is continuous at {pc,S) if 
6 > Pc- e). 

Finally, we consider the case where p > Pc- If 9{p, 5) = 0, continuity at 
{p,6) follows from part (b) of Proposition 3.3. So suppose 9{p,6) > 0. Let 
a as in Theorem 5.1. By that theorem there is an n > N{p) with 

fp,s{^,n) > 1 - a. 

Hence, by Lemma 3.1 there is an open neighborhood W of (p, 5) such that 

fp>,S'{S,n)>l-a, (38) 

for all {p',6') € W. Since n > N{p) and N{.) is decreasing, it follows 
from (38) and Theorem 5.1 that 9{-, •) > on S, where S is the set of all 
{p', S') G W with p' > p. From this and Lemma 2.2 we conclude that 9{-, •) 
is also strictly positive on the set 

U := {{p',S') : <p and p' + {l-p')S' > r + {l-r)(3 for some {r,/3) G S}. 

It is easy to see that S UU contains an open neighborhood of {p,S). Now 
it follows from part (a) of Proposition 3.3 that 9{-) is continuous at ip,5). 
This completes the proof of the main theorem. □ 

Remark 6.1. A crucial role in the proof is the finite- size criterion, Theorem 
5.1. That theorem has been formulated for p > Pc. When p = Pc (or < Pc) 

the sdp model is an ordinary percolation model, for which a similar criterion 
is known. In fact, for ordinary percolation we do not have the dependency 
problems which led to the introduction of N . Consequently, for p = Pc 
we can take N = 1. But, on the other hand, if we let p tend to Pc from 
above, the upper bound on N{p) obtained from our arguments in Section 
5 tends to oo. And that, in turn, comes from the fact that our bound on 
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dependencies, Lemma 2. 3, is in terms of path probabilities for an ordinary 
percolation model ( on the matching lattice, with parameter 1 — p) which is 
subcritical but approaches criticality (which makes these bounds worse and 
worse) as p approaches Pc from above. This is essentially why the proof of 
Theorem 1.2 does not work atp^. Of course, if it would work, the conjecture 
referred to in Section 1 would be false. We hope that attempts to stretch the 
arguments in our paper as far as possible will substantially increase insight 
in the conjecture and help to obtain a solution. 
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